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Abstract. We study the dynamics of a complex open quantum many-body system. 
The coupling to external degrees of freedom can be viewed as a coupling to a radiation 
field, to continuum states or to a measuring apparatus. This perturbation is treated 
in terms of an effective non-Hermitian Hamiltonian. The influence of such coupling on 
the properties of the many-body dynamics is discussed, with emphasis on new effects 
related to dynamical segregation of fast and slow decays and the phase transition to 
Dicke superradiance. Relations to quantum optics, continuum shell model, theory 
of measurement, quantum chaos, percolation theory, and to quantum reactions are 
stressed. 
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1. Introduction 

The interplay of complex intrinsic dynamics with radiation, coupling to the continuum 
or to any other external influence, such as measurement, is an important subject for 
quantum optics, nanoscale devices, theory of decoherence and quantum complexity, 
future quantum computing, and physics of nuclei far from stability, to name just few 
branches of modern physics, seemingly far remote from each other. After all, the 
complete isolation of a "system" and the "external world" is never possible. The chaotic 
and regular aspects of many-body physics, as well as openness and decay properties, are 
deeply rooted in valuable parts of Wigner's legacy. In this paper, commemorating 
Wigner's contributions, we would like to stop at the crossroad of ideas related to 
the complexity of many-body physics and the Weisskopf-Wigner damping theory. We 
consider a mesoscopic many-body system with internal dynamics governed by the mean 
field and residual two-body interactions. We assume that the system is perturbed by 
non-Hermitian terms in the Hamiltonian. These terms, in general, model an effective 
coupling to the external degrees of freedom such as the continuum of scattering states, 
radiation fields or measurements. 

The subject of this work is cross-disciplinary and the approach implemented here 
is widely used in atomic, molecular and nuclear physics, see our recent work pQ and 
references therein. For a broader audience we structure the paper and conduct our 
discussion in a pedagogical manner. We start our work by presenting in Sec. El a well 
known description of the internal dynamics in an open many-body quantum system with 
the aid of the non-Hermitian and energy- dependent Hamiltonian. A simple example of 
two interacting spins in the presence of decay is discussed in Sec. H3 Although this 
example has no complexity associated with the many-body dynamics, the important 
physical features related to decay, as well as mathematical aspects, are already there. 
This case serves as an introduction to our main discussion presented in Sec. HJand based 
on the consideration of a system of interacting fermions coupled to a decay channel. 
Concluding remarks are given in Sec. El 

2. Effective Hamiltonian of an open system 

The technique using a non-Hermitian and energy-dependent effective Hamiltonian is 
rather common and dates back to the Weisskopf-Wigner damping theory |2j, works of 
Rice (3] and Fano |3J, and Feshbach projection formalism j^]. It is assumed that the 
many-body states in the system can be separated into two classes, the internal (intrinsic) 
states and external ones. The set of internal states |a) is generally understood to be 
composed of the states where all particles occupy bound mean-field orbitals. The fact 
that the system is open requires one to include into consideration a set of external 
many-body states. These states \c; E) are viewed as decay or reaction states and contain 
particles in the continuum or radiation quanta; we label them with a continuous variable 
of energy and a set of other quantum numbers characteristic for a channel c. Discussion 
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of the reactions is not our primary goal in this work, nevertheless we should emphasize 
Wigner's significant contribution to this subject 00 IE]- 

Although there exist many different formulations of this technique, here we start 
from the Schrodinger equation 

H\V) = E\V), (1) 

where the full wave function is a superposition of internal and external states, 

\^} = J2C a \a) + J2 f dE X c E \c;E). (2) 

a c 

Eliminating the external states we come to the equation for the coefficients Cp that 
describe the internal part of the wave function, 



C„ = 0. (3) 



This equation looks like a traditional eigenvalue problem with the effective Hamiltonian 
matrix in the intrinsic space defined as 

(a\H{E)\(3) = (a\H\0) +£) f dE ' ^^^l , (4) 

c J 

where we introduced A c a (E) = (a\H\c;E), the coupling amplitudes between internal 
and external states. 

The first term in Eq. (j3J) is the standard Hamiltonian matrix that would describe 
our system if it were decoupled from the external world. For definiteness, we assume 
that this part of the Hamiltonian is 

H = ^ei4ai + ^12;34aI4 a 3 a 4 , (5) 

1 1234 

with single-particle energies e\ and residual two-body interactions with antisymmetrized 
matrix elements Pi2;34 • The indexes 1, . . . , contain a complete set of single-particle 
quantum numbers, and the total many-body basis state a can be constructed, in the 
case of Fermi-statistics, as a Slater determinant of these single-particle states. 

The integral in Eq. (|IJ) can be further decomposed into its Hermitian part (principal 
value) and the remaining non-Hermitian part, 



, Al{E')A c JE 



E-E' 



' E _% -™ E wmm ■ ( 6 ) 

c (open) 

Thus, the effective Hamiltonian (J7j) can be written as 

H(E) = H + A(E)- t -W(E), (7) 



n h a 



w 



Figure 1. A schematic representation of the effective Hamiltonian of an open system. 
The Hermitian shell model part H is shown by the two terms following Eq. J5|, single- 
particle energies and two-body residual interactions. The term A includes all channels 
c while the term {—i/2)W is restricted to the channels open at a given energy. 



which contains an energy-dependent self-energy term due to the "off-shell" processes of 
virtual excitation into all, open and closed at given energy, continuum channels, 



The imaginary part W has a factorized form, and originates from the real "on-shell" 
processes through the continuum channels that are open at given energy. The schematic 
diagrams representing various terms in (JJJ) are shown in Fig. ^ 

For our current discussion we will ignore the energy dependence of the Hamiltonian 
(J7J). In general, the energy dependence of A is smooth and frequently this Hermitian 
term can be incorporated into the intrinsic Hamiltonian if as a renormalization. The 
energy dependence of the non-Hermitian part W can be crucial, in particular near 
channel thresholds, see for example Wigner's work jH]. In loosely bound systems this 
term determines the exact binding energy p. Here, however, we are going to discuss 
the new effects and especially the phase transition generated by the presence of non- 
Hermitian terms in the Hamiltonian, and for this purpose the amplitudes A c a can be 
taken as energy-independent parameters, as it would take place in a physical situation 
with remote thresholds. We simply assume that the Hamiltonian TC = H — iW/2 is 
known, the ant i- Hermitian part has a factorized structure (J9jl . and this will allow us to 
study the resulting dynamics. 

3. Two spin-system 

We start with a very simple example similar to the system discussed by Dicke jH] in 
relation to superradiance. Consider two interacting distinguishable spin-1/2 molecules, 
Si = S2 = s = 1/2 , with the spin-spin interaction 





and explicitly non-Hermitian and energy dependent term 




(9) 



c (open) 



H° — a Si ■ s 2 . 



(10) 
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This Hamiltonian is diagonal in the basis of the total spin S = Si+s 2 and its z-projection, 
and the energies of the states, 

E s = ^[S(S+l)-2s(s + l)) , (11) 

are E = — 3a/4 for the singlet, S — 0, state and E\ = a/4 for the triplet, 5 = 1, state. 

Now let us place the system in the magnetic field that produces two effects. First, 
it leads to the Zeeman splitting that can be described by the additional term in the 
Hamiltonian 

H B = es\ + es z = eS z . (12) 

For simplicity we assume that the field acts in the same way on both molecules, therefore 
the additional term (fT2"|) still commutes with the operator S 2 . Secondly, we assume that 
this two-spin system in the magnetic field becomes open; in the presence of the field 
the first molecule in its excited polarized state can dissociate, the phenomenon similar 
to the Feshbach resonance [ID] . This means that the molecule in sf — 1/2 state decays 
exponentially with half-life T1/2 which can be described by the width 7 = In 2/7^. 
In terms of the general formalism of Sec. El here we have one open decay channel with 
coupling strength A and 7 = A 2 . In principle the decay rate of a molecule should 
depend on the magnitude of the Zeeman splitting and thus on the field strength. For 
the interacting system the decay clearly should depend on total energy. As we stated 
earlier, we ignore this dependence in our examples. The situation can be modelled by 
an extra non-Hermitian term in the Hamiltonian 

W = -i 1 -{s z 1 + s), (13) 

so that the first molecule in the state with s\ = 1/2 would have decay width 7 and 
this width would be zero in the state s\ = —1/2. As a result, in the magnetic field the 
effective Hamiltonian for the two-spin system becomes 

W = asi-s 2 + esi + esf-»^(si + s). (14) 

The breakup of the first molecule distorts the symmetry of the system, S 2 no longer 
commutes with 7i, and the singlet and triplet states become mixed. The situation 
here is interesting in the sense that the open decay channel introduces an incompatible 
symmetry: the interaction Hamiltonian ([TTj) and the magnetic field part ()12|) are both 
diagonal in the basis of the total spin (5* , S z , si , s 2 ), whereas the decay part, as that 
of an intrinsic elementary process of the molecule 1, is diagonal in the decoupled spin 
states of molecules (si , s\ ; s 2 , s 2 ), the so-called m-scheme basis. The axial symmetry 
with respect to the magnetic field holds for the full Hamiltonian (jl4j) requiring the 
conservation of S z . 

The full effective Hamiltonian matrix of the S z = block in the m-scheme basis 
takes the form 

n = -i + U ;; 1. do) 
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Figure 2. Evolution of two states £ = E — iT/2 in the complex plane as a function 
of the increasing parameter 7. 

and the diagonalization of this matrix determines the quasistationary eigenstates of the 
system and their complex energies, 

The trajectories of the complex energies S± as functions of the increasing decay 
parameter 7 are shown in Fig. El 

The dynamics of the two-level effective Hamiltonians have been extensively studied 
in the past, see for example [Hl[T2l[T31IIlj. In fact the discussion was initiated by Wigner 
and von Neumann ^5], who demonstrated that mixing by an off-diagonal Hermitian 
interaction leads to the level repulsion. The behavior of complex eigenstates as a function 
of real and complex parts of the effective Hamiltonian can be summarized as follows. It 
is well known that the increase of the Hermitian part of the interaction leads to repulsion 
of real energies E and attraction of widths T since the decay probability is shared by 
the two mixed states. It is less well known that the non-Hermitian part has an opposite 
effect, it increases the difference between widths and moves real parts E closer together 
[TTJ H2 IH] ; these effects are schematically demonstrated in Fig. |3] 12 . 

Our example reveals the generic features, consider Eq. (jT5j) and Fig. El In the 
region of weak decay, i.e. slow process of molecular disintegration, < 7 < 2a, the two 
eigenstates £ + and £_ have equal widths, T + = T_ = 7/2, and the real energies E + and 
E- are attracted to each other as 7 grows. At the critical point 7 = 2a the levels cross 
and the "phase transition" occurs as a sharp change in the behavior of the two complex 
energies. The real parts become degenerate, E + = E_ = —a/4, while the widths 
separate, T± = j/2 ± a/ (t/2) 2 — a 2 so that in the limit of strong decay, 7 > a, we have 
a single "superradiant" state with the width r + — > 7 accumulating the entire available 
decay probability, and one long-lived state with r_ — ► 0. The f actor ized structure of 
matrix W (JHJ), that results in zero eigenvalues in the limit when this part dominates the 
Hamiltonian, is a key feature responsible for this behavior. As can be seen from the 
derivation of Sec. 121 this structure is related to the analytical structure of the problem 
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Figure 3. Generic behavior of two resonances E± and their widths T± as a result of 
increase of the Hcrmitian part of mixing interaction V and in response to increase of 
non-Hcrmitian mixing W. 



and unitarity requirements. The formation of fast and slow decay modes, which can be 
viewed as an ultimate width repulsion, is a typical effect in systems with overlapping 
resonances jTHJ [T2j . In particular, a detailed discussion for two-level examples can be 
found in ITH| ITU H%] and references therein. 

The coupling to decay, as we see, has important observable consequences for both 
the intrinsic state of the system and the external scattering picture. For example, for 
a singlet ground state (a > 0), the magnetic field does not affect the system in the 
absence of molecular breakup, because Hb is identically zero in the singlet state. In the 
presence of the field, the disintegration of the first molecule reorients the eigenvector 
of the system admixing a triplet state. Thus, the system becomes susceptible to an 
excitation by the magnetic pulse. 

Looking at the situation from a different angle, namely from the "exterior", it is 
clear that the decay of a first molecule alone is different from the decay properties 
when the second molecule is present. In the first paper on superradiation jU] Dicke 
argued that coupling to the radiation field can lead to strong modification of the 
intrinsic properties of the system. In context of nuclear reaction theory, a significant 
contribution to this subject was made by Wigner and his collaborators who discussed 
the nature of resonances in the nucleon-nucleus interactions, sum rules related to the 
widths, distribution of the widths among the states and general features of the scattering 
cross section [HI HH] , that in context of our example can be related to the properties 
of the matrix (|15j). 

4. A model of an open many-body system 

The transition to superradiance in a many-body model with decay and random internal 
dynamics is our second example. We consider an equidistant set of Q fermionic single- 
particle levels ( "orbitals" ) and assume that the upper level, that we label as u, belongs to 
the continuum, i.e., in the approximation of non- interacting particles, the single-particle 
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Figure 4. Schematic drawing of the single-particle states in the model of Sec. 0| The 
upper state belongs to the continuum. 



wave function of this orbital decreases exponentially with time, while the remaining 
orbitals are particle-stable, see Fig. EJ This situation again corresponds to a single 
decay channel where any many-body state with an occupied single-particle orbital v is 
coupled to the continuum with the corresponding amplitude A. Equivalently the full 
non-Hermitian effective Hamiltonian TC can be constructed from (jHJ) under assumption 
that one of the single-particle energies is complex, 

e* = e„-|7, 7 = ^ 2 ; (17) 

it contains single-particle energies with the scale determined by their spacing Ae , the 
non-Hermitian part W, the scale of which is determined by 7, and the Hermitian mixing 
V, Eq. (jSJ). We do not impose any symmetry here; the effects of Kramers degeneracy 
and rotational invariance have been discussed in pQ. For our example we select the 
parameters of the Hermitian residual two-body interaction V, eq. (JHJ at random from 
the Gaussian distribution centered at zero. The average absolute magnitude of these 
matrix elements \fv^ sets the energy scale for the model. In what follows we assume 
\J~V^ = 1 and allow Ae and 7 to vary. 



4-1. Evolution of many-body states 

A system of N fermions in Q single-particle states has 



M 



n\(q-n)\ v ' 

many-body states. To be specific, we consider iV = 4 particles in Q = 8 single-particle 
orbitals. The M = 70 eigenvalues Sj = Ej — (i/2)Tj of the Hamiltonian are in the lower 
part of the complex plane, and their motion as a function of the increasing parameter 
7 is shown in Fig. here Ae = 1. This figure is a generalization of the situation in 
Fig. discussed in the previous section El As stated earlier, it is assumed that 7i is 
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independent of energy E . Therefore the trace conservation leads to a helpful sum rule 
for the imaginary part of the Hamiltonian, 

Vr,( 7 ) = -2Im(TrH). (19) 



j 



This sum rule is a direct analog to the sum rules discussed by Wigner and collaborators 
in the context of neutron resonances [0]. 

Several dynamical limits can be identified. In the "shell-model" limit of 7 = 
and TC — H, standard diagonalization results in orthonormal eigenstates ^.m. with real 
eigenvalues E sra . In Fig. the values E sra define the starting points on the real axis. 
As soon as the single-particle width 7 7^ 0, the decay channels open. The random 
residual interaction V leads to a chaotic mixing of configurations, and all many-body 
states acquire widths. This excludes special cases when configurations with the occupied 
decaying state are blocked, for example by selection rules. 

For small 7, the term W in the Hamiltonian can be treated perturbatively. In the 
lowest order this results in 

£ = Km. ~ \ (%.m\W\V 8 . m ) . (20) 

Then the decay width of a many-body state is determined by the occupancy n v (^>) = 
(^ s . m \ala u \ty s , m ) of a decaying orbital v calculated for a particular shell model state 

* B .m., 

iXtt) = (* B .m|W|*B.m) = 7M*s.m.) • (21) 

Only in this limit it is possible to establish the "natural" result (121 J) that is widely used 
in extracting spectroscopic factors of intrinsic states from various reaction cross sections. 

The opposite limit of strong 7 is less trivial. It is clear from Fig. 1 that all many- 
body states are divided into two groups. In the first group the complex energies have 
a rapidly increasing width, whereas the states in the second group become long-lived 
with their widths approaching zero. The dominance of the external coupling W is the 
origin of this dynamical separation. The single-particle energies and mixing V are of 
minor importance in this limit, and the states can be classified by their relation to the 
decay. The matrix W is highly degenerate; because of its factorized structure, Eq. 0, 
it has only 

\r _ ("-!)! (22) 
Mr ~ (N-1)W-N)\ (22) 

nonzero eigenvalues. The corresponding eigenstates of W are the configurations that 

have the decaying upper orbital occupied and thus acquire in the strong decay limit 

the maximum width equal to 7 so that the trace (J19)) of the imaginary part of the 

Hamiltonian is exhausted by jNr- 111 our case, the total dimension of the many-body 

space is M = 70, and Ar = 35. The remaining 35 states correspond to the empty upper 

orbital, and correspondingly to the zero eigenvalues of W, so that decay of these states 

is strongly hindered. 
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Figure 5. Complex plane trajectories of all 70 many-body states in the system of 
4 particles in 8 equidistant single-particle levels with the spacing of Ae = 0.5 energy 
units as a function of the decay width 7 of the upper single-particle state increasing 
downwards. The two-body interaction matrix elements are taken randomly from the 
Gaussian distribution with zero mean and variance of one energy unit. 



Thus, in the limit of strong coupling to the continuum, a dynamical segregation of 
short-lived and long-lived (compound) intrinsic states occurs. This is a generalization 
of the effect observed in a simple two-level system of Sec. H3 This phenomenon was 
discovered in numerical simulations (201 EI], see also [221123], was studied in molecular 
and atomic physics [21] |2HJ I2E] and observed in experiments with microwave cavities 
|27j . The fact that the underlying physics is analogous to the Dicke coherent state j§] 
was explained in Refs. [TT)] I17j. The coupling of states via a common decay channel 
(width collectivization) plays a role similar to that of the common radiation field of 
spontaneously radiating atoms for the Dicke superradiant coherent mode experimentally 
observed in quantum optics [23 EH] ■ 

A related effect considered in the quantum theory of measurements has possible 
ramifications for quantum computers. An external influence on the system brought by 
a procedure of measurement, the Zeno effect [30], can be modelled similarly, with the 
aid of non-Hermitian terms in the Hamiltonian. Recently this topic attracted a lot of 
attention due to the appearance of remarkable experimental data [HUJ EH] 02] ■ in this 
context it was demonstrated by Facchi and Pascazio [HB] that an analogous effective 
subspace separation plays a vital role. We need to emphasize that the dominating 
term in the Hamiltonian, W, because of its degeneracy, only effectively decouples the 
subspaces but it does not fully determine the dynamics. The single-particle energies 
and residual Hermitian interactions determine the motion within the intrinsic subspace 
that still can be chaotic. 

The situation at intermediate values of the decay strength, as seen from Fig. is 
complicated. In this region all three terms of the effective Hamiltonian are important, 
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and their interplay is responsible for chaotic transitional dynamics. As 7 increases, 
the perturbative approach breaks down. However, at this point the limit of the full 
subspace decoupling is not yet reached. In this regime we are faced with a peculiar 
phase crossover, where all components of the dynamics are strongly mixed and thus 
both energy centroids Ej and widths Tj become affected in a non-trivial way by the 
presence of the non-Hermitian component. A similar interplay between single-particle 
energies and the residual pairing interaction is known to result in chaotic dynamics of 
the normal-to-superconducting phase transition in mesoscopic systems [31] which, due 
to the finite nature of the system, occurs as an extended crossover. 

In his work related to nucleon- nucleus interaction Wigner emphasized the 
importance of interplay between the state of intrinsic configuration (determined by 
Hermitian interaction V) and doorway states (the eigenstates of W). He pointed out 
that three types of behavior can be distinguished by the distribution of the widths among 
many-body states. In the limit of V — > 0, the "independent particle model," the entire 
width is absorbed by the set of states where the decaying single-particle orbital is fully 
occupied. This is an analog to the superradiant regime. In the intermediate regime, 
due to the interaction V, the width becomes fragmented and spreads to many states. 
Finally, in the limit of large V the chaotic intrinsic motion leads to the approximately 
uniform distribution of the widths, the "uniform model," which is analogous here to the 
shell-model limit 7^0. 



4-2. Effective occupation numbers 

We can extend the definition of the single-particle occupancy n u of the decaying orbital 
to any value of 7, instead of the 7 — > limit discussed previously, eq. PTjl. The effective 
occupation number for the quasistationary state j can be defined as 

nAr<7) = 9 IM. (23) 

As follows from Eq. (JTHJ), the sum Yjj n vU'i 7) is independent of 7 . The numbers (J23J), 
properly normalized according to Eq. (JEJ), show how additional decay amplitudes 
induced by an infinitesimal increase of 7 are distributed among the widths of the 
eigenstates. Despite the resemblance to occupation numbers, these n v can be negative 
and do not have to be confined to the interval between and 1. There exists also 
another possibility, namely, to introduce spectroscopic factors directly as 1^(7) /7 . They 
are indeed between and 1, but being simply proportional to T they do not reveal any 
new information related to the dynamics. 

In Fig. El the evolution of the effective occupation numbers n v for the lowest 13 
states tyj is shown. The evolution starts with occupancies found for the shell model with 
residual interaction but no decay and follows a complicated path to the final segregated 
values of and 1, corresponding to quasistable and superradiant states, respectively. 
The resulting crossover transition of the entire spectrum is gradual. It starts at the 
point where 7 exceeds the level spacing on the real axis. Here the resonances start 
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Figure 6. The evolution of the effective occupation numbers n v , Eq. (|23|l . as a 
function of 7 for the 13 lowest (selected at 7 = 0) many-body states. 

to overlap. Although there are states that converge to the limiting value relatively 
fast, the majority of states undergo a complicated and peculiar evolution driven by 
the competition between the intrinsic and external interactions. Before reaching the 
regime of developed segregation, the evolution may involve peaks of occupation numbers 
for specific eigenstates. The segregation occurs at the values of width covering the 
macroscopic portion of the many-body spectrum, so that a balance between the external 
coupling represented by the imaginary part and the internal motion given by Hermitian 
part is reached. An analogous balance between two-level atoms (radiators) and radiation 
field is required for the Dicke superradiance. The nature of the phase transition to 
superradiance in the Dicke model was studied in Refs. jH3 EE] ■ 

The definition of the effective occupation numbers in the form (|2*Hj) provides a 
special insight into the dynamics of a phase transition. This definition is similar in 
spirit to that of correlational entropy |HZ| that shows how fast the states of the system 
readjust as a result of changes in the environment. Here 7 plays the role of the external 
parameter. For regular dynamics the behavior of this entropy is smooth. However, 
any phase transition, which may be hindered in a mesoscopic system, leads to a rapid 
restructuring of the state and thus results in a peak in the correlational entropy jHH EEj • 
Similarly, in our example, in the transitional region of 7, the peaks of Thy cLS cL function 
of 7 are clearly seen for most of the states ^j. For some states in Fig. 2 the peaks 
can be sharp, an indication that a particular many-body state undergoes a sudden 
restructuring. It is worth noting that a practically identical discussion is appropriate 
for the two-spin example considered in Sec. El There the quantity n v becomes the 
effective probability of finding the molecule 1 in the radiating state with s z = 1/2. For 
the two mixed states with S z = it can be easily calculated using the definition 
and the expression for the complex energy (jlfijl . Even in the limited example of two 
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Figure 7. Fraction of segregated states as a function of 7 . The single-particle level 
spacing is varied from the degenerate case, Ae = 0, to Ae = 100, a point where the 
residual interaction V can be completely ignored. 



spins the "phase transition" (the point of level crossing) appears as a special bifurcation 
point. 

4-3. Global transition to superradiance 

In Fig. 0we demonstrate the global properties of the transitional regime from another 
perspective by plotting the segregation fraction £(7) defined as a fraction of many-body 
states for which the generalized occupancy n u appears close to the limiting values of 
or 1, 

£( 7 ) = — ^ J e - n ' (i;7)/{2<T2) + e- [1 ~ nlij]l)]2/{2a2) ~ . (24) 
3 

Here we choose a = 0.1. The role of the relative magnitudes of the single-particle 
spacing Ae and two-body matrix elements for the phase transition can be inferred 
from this figure. Here a randomly chosen set of two-body matrix elements V with an 
average magnitude setting the unit of energy is kept constant, while the single-particle 
level spacing Ae is varied. In the limit of infinitely large spacing, the mixing part of 
the Hamiltonian, V, becomes negligible, formally resulting in the disappearance of the 
intermediate transitional region. Nevertheless, even for Ae fifty times larger than the 
average V, the transitional region is still well defined. 

The transition to superradiance, being gradual for the entire system, can be very 
drastic for a particular many-body state. An interesting analogy can be found between 
the transition to superradiance and percolation theory. The set of original intrinsic 
states can be considered as a lattice with cites determined by energies E s _ m _. Without 
decay the system is stationary and rests at a given initial eigenstate. The coupling to 
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the continuum opens the possibility for the system to hop from state to state until it 
gets into a doorway state and emits a particle. The probability of such hopping j — > j' 
according to the perturbation theory is determined by the matrix element (j'lW^lj') (this 
term in our examples is proportional to 7) and inversely proportional to the distance 
l-E's.m.j — E 8 . m .j,\. The lattice is irregular because the cites of the lattice correspond 
to the eigenvalues of the Hamiltonian H and the distribution of distances between 
the cites is given by Wigner random matrix theory [55] . As a result, some areas can 
experience phase transitions and coherently orient toward a doorway state (which is 
not an eigenstate) earlier than the rest of the system. This is the reason for the sharp 
changes observed in individual states, see Fig. [H] and related discussion. As follows from 
Fig. [7[ a significant strength of the continuum coupling is required for the full system 
to reorient and to create a coherent superradiant flow into doorway configurations. 

5. Conclusions 

We presented an introduction to interesting phenomena that appear in complex systems 
as a result of their coupling to the environment. We discussed simple examples showing 
that the strong coupling significantly reorganizes the motion, dynamically separating the 
states by their relation to the external world. We investigated the properties of a many- 
body phase transition from the normal shell model (or Fermi-liquid) behavior of a closed 
many-body system to the superradiant regime with the segregation of rapidly decaying 
states from long-lived resonances. The chaotic character of the intrinsic interaction 
does not prevent the system from this generic transition. The presence of very long- 
lived states, that are almost decoupled from the external world, maybe important for 
the problem of quantum computing. 

We deliberately related our examples to different branches of physics and provided 
many references trying to demonstrate that these phenomena can be relevant to almost 
any case when the coupling of a system to the outside world cannot be neglected. 
This topic links together many of the ideas inherited from Wigner. They include 
the complexity of many-body states with the characteristic level repulsion [15] ; the 
theory of quantum damping and decay [2|; the role of symmetry and incompatibility 
between symmetries represented here by the competition between internal symmetry 
and symmetry imposed by the decay processes; theory of resonances, resonance sum 
rules, fragmentation of the widths in reactions [HI El El and the near-threshold width 
behavior [Hj; random matrix theory [55] . 
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